Nearly Perfect Single-Channel Conduction in Disordered Armchair Nanoribbons 
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The low-energy spectrum of graphene nanoribbons with armchair edges (armchair nanoribbons) 
is described as the superposition of two non-equivalent Dirac points of graphene. In spite of the lack 
of well-separated two valley structures, the single-channel transport subjected to long-ranged impu- 
rities is nearly perfectly conducting, where the backward scattering matrix elements in the lowest 
order vanish as a manifestation of internal phase structures of the wavefunction. For multi-channel 
energy regime, however, the conventional exponential decay of the averaged conductance occurs. 
Since the inter- valley scattering is not completely absent, armchair nanoribbons can be classified 
into orthogonal universality class irrespective of the range of impurities. The nearly perfect single- 
channel conduction dominates the low-energy electronic transport in rather narrow nanorribbons. 

PACS numbers: 72.10.-d,72.15.Rn,73.20.At,73.20.Fz,73.23.-b 
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I. INTRODUCTION 

Graphene is the first true two-dimensional (2D) ma- 
terial 1 . Due to the honeycomb lattice structure of sp 2 
carbon, the n electronic states near the Fermi energy 
behave as the massless Dirac fermions. This leads to 
many nontrivial properties of graphene such as half- 
integer quantum Hall effect 2 . The valence and conduc- 
tion bands touch conically at two non-equivalent Dirac 
points, called K+ and K.— points, which possess oppo- 
site chirality 3 -. In graphene, the presence of edges can 
have strong implications for the electronic band struc- 
ture of 7T electrons^££. Graphene nanoribbons (GNR) 
with zigzag edges are known to have partial flat bands 
near the Fermi energy due to the edge localized states. 
The electronic structures of nanoribbons with armchair 
edges crucially depend on the ribbon width 4,5,6 - 7 . Re- 
cent rapid progress of experiments confirmed the edge- 
dependent electronic states of graphene using scanning 
tunneling microscope^, and also succeeded in creating 
GNR using lithographic^ or chemical techniques 1 ^. 

GNR displays unusual electronic transport properties, 
in apparent conflict with the common belief that ID sys- 
tems are generally subject to Anderson localization. In- 
deed it was demonstrated that nanoribbons with zigzag 
edges (zigzag nanoribbons) with long-ranged impurities 
possess one perfectly conducting channel (PCC), i.e. the 
absence of Anderson localizatio n 12 : 13 . Since in zigzag 
nanoribbons the propagating modes in each valley con- 
tain a single chiral mode originating from edge states, a 
single PCC emerges, associated with such a chiral mode, 
if the impurity scattering does not connect the two val- 
leys, i.e. for long-ranged impurities (LRI). 

In this paper, we show that the single-channel trans- 
port in the disordered armchair nanoribbons subjected to 
the long-ranged impurities is nearly perfectly conducting 
in spite of the lack of well-separated two valley struc- 
tures. The origin of the nearly perfect conduction is the 
cancellation of the backward scattering matrix elements 
in the lowest order due to the manifestation of internal 
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FIG. 1: (a) Structure of graphene armchair nanoribbon. The 
area with the length L represents the disordered region with 
randomly distributed impurities, (b) Energy dispersion of 
armchair ribbon with N = 14. The energy range for single- 
channel transport is described by A s . (c) 1st Brillouin Zone 
of graphene. 



phase structures of the wavefunction. For multi-channel 
energy regime, however, the conventional exponential de- 
cay of the averaged conductance occurs. Since the inter- 
valley scattering is not completely absent, the disordered 
armchair nanoribbons can be classified into orthogonal 
class. The nearly perfectly conducting effect dominates 
the low-energy electronic transport properties in rather 
narrow nanorribbons. 

The paper is organized as follows: In Sec. [Til the tight- 
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FIG. 2: (a) Average conductance (g) as a function of the ribbon length L in the presence of long-ranged impurities (LRI) 
for several different Fermi energies E. Conductance is almost unaffected by impurities for single-channel transport (E = 
0.1,0.2 and 0.3) while it shows a conventional exponential decay for multi-channel transport (E > 0.4). Here, TV = 14, 
JT-imp. = 0.1 and d/a = 1.5. Ensemble average is taken over 10 4 samples, (b) The Fermi energy dependence of (g) for LRI. (c) 
The same as (a) for short-ranged impurities (SRI). Here, TV = 14, nj mp . = 0.1 and d/a = 0.1. 



binding model used in our numerical simulation is ex- 
plained. We also briefly review the electronic states of the 
low-energy single channel mode in armchair nanoribbons 
by k ■ p scheme. In Sec. lIIIl we present the numerical 
results indicating the nearly perfect single-channel con- 
duction. This property is then explained by T-matrix 
analysis. Symmetry consideration is also given in this 
section. Finally we summarize our work in Sec. lIVI 



II. 



ELECTRONIC STATES OF ARMCHAIR 
NANORIBBONS 



Gaussian form of range d, 

V l = V(r i )= uex P 
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where the strength u is uniformly distributed within the 
range |u| < um- Here um satisfies the normalization 
condition: um 



(funspacc) exp( _ r 2 /d2)/( ^3 /2) = ^ 



this work, we set m, 



0.1, wq = 1.0 and d/a — 1.5 for 



LRI and d/a = 0.1 for short-ranged impurities (SRI). 



A. Tight-binding model 



B. Low-energy single channel mode 



We describe the electronic states of graphene nanorib- 
bons with armchair edges by the tight-binding model 



(i) 



where c$(cj) denotes the creation (annihilation) opera- 
tor of an 7r-electron on the site i neglecting the spin de- 
gree of freedom. 7^ = — 1 if i and j are nearest neigh- 
bors, and otherwise. In the following we will also ap- 
ply magnetic fields perpendicular to the graphite plane 
which are incorporated via the Peierls phase: 7^ — > 



lij exp 



i2ir(e/ch)fi dl ■ A 



where A is the vector po- 



tential. The second term in Eq. JT]) represents the impu- 
rity potential, V% = V(r{), at position r». 

In Fig.[T]Ja), the schematic figure of armchair ribbons is 
depicted. The ribbon width TV is defined by the number 
of zigzag kinks in the transverse direction. The armchair 
ribbon can be metallic if TV — 3m— 1 (to: integer number) 
as shown in Fig.[ljb), otherwise semiconducting. The 
disordered sample region with the length L is attached 
to two reservoirs via semi-infinite ideal regions. 

We assume that impurities are randomly distributed 
with density nj mp .. Each impurity potential has the 



Here we briefly review the relation between the low- 
energy electronic states of armchair nanoribbons and the 
Dirac spectrum of graphene. The electronic states near 
the Dirac point in graphene can be described by the mass- 
less Dirac Hamiltonian 



Ho = 7 k x {<J x ® r°) - k y (a y (g) t z ) 



(3) 



acting on the 4-component pseudospinor envelope func- 
tions F(r) - [F+(r), F+(r), Fj(r), (r)] , which 
characterize the wave functions on the two crystalline 
sublattices (A and B) for the two non-equivalent Dirac 
points (valleys) shown in Fig.[]Jc).The correspond- 
ing wave vector for the K+ point is If = (27r/a)(2/3, 0), 
and that for the K_ point is —K. We have defined 
the amplitude of wavefunction at Ra and that at Rb 
as iPa{Ra) = e iKRA F+(R A ) + e- iKRA FX(R A ) and 
iPb{Rb) = e iK RB F+{R B ) - e- iK RB Fs(R B ), respec- 
tively. Here, Ra (Rb) is the coordinate of an arbitrary 
A(B) sublattice site. Here 7 is the band parameter, 
k x (ky) are wave number operators, and t° is the 2x2 
identity matrix. Pauli matrices a x ' y ' z act on the sublat- 
tice space (A, 73), while t x ' v,z on the valley space (K±). 
The boundary condition for armchair nanoribbons^ can 
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be written as 



[F+(x,y) + FX(x,y)]\ x=0 , w =0, (4) 
[F+(x,y)-F B {x, V )]\ x=0 . w ={). (5) 

Since this boundary condition projects and _KT_ 

states into T point in the first Brillouin Zone as seen in 
Fig.fljc), the low-energy states for armchair nanoribbons 
are the superposition of K+ and K_ states. If the ribbon 
width W satisfies the condition of W = (3/2)(N w + l)a 
with N w = 0, 1, 2, ... , the system becomes metallic with 
the linear spectrum. The corresponding energy is given 
by 



(6) 



where k t 



0,±1,±2,... and s = ±. 



3(N w + l)a ' 

The n = mode is the lowest linear subband for metallic 
armchair ribbons. The energy gap (A s ) to first parabolic 
subband of n = 1 is given as 



Da, 



(7) 



which is inversely proportional to ribbon width. It should 
be noted that small energy gap can be acquired due to the 
Peierls distortion for half-filling at low temperature o 7 ' 15 , 
but such effect is not relevant for single-channel transport 
in the doped energy regime. 
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FIG. 3: Inverse of localization length £~ as a function of the 
magnetic flux <j> through a hexagon ring measured in units of 
ch/e, in the presence of (a) LRI and (b) SRI for different 
Fermi energies E. In (a), £ _1 for E = 0.4 is omitted since its 
value is much larger than others. 



the dimensionless conductance by using the Landauer 
formula, g{E) — Tr(tH). Here the transmission matrix 
t(E) for disordered system is calculated by using the re- 
cursive Green function method^. 

Figure [2a) shows the averaged conductance (g) as a 
function of the ribbon length L in the presence of LRI 
for several different Fermi energies E. As we can clearly 
see, the averaged conductance subjected to LRI in the 
single-channel transport (E = 0.1,0.2 and 0.3) is nearly 
equal to one even in the long wire regime. This result is 
contrary to our expectation that electrons are scattered 
even by LRI, since wavefunctions at K+ and _K_ points 
are mixed in armchair ribbons. For multi-channel trans- 
port (E > 0.4), the conductance shows a conventional 
decay. The robustness of single-channel transport can be 
clearly viewed from the Fermi energy dependence of con- 
ductance for several different ribbon lengths L as shown 
in Fig.^b). It should be noted that the energy depen- 
dence in the vicinity of E = is quite different from that 
in zigzag nanoribbons. The conductance decays rapidly 
due to the finite ribbon width effect in zigzag ribbons 12 
while the conductance around E = remains unity in 
armchair ribbons (Fig.[^b)). 

Now let us see the effect of short-ranged impurities 
(SRI). Figure [DJc) shows the average conductance (g) as 
a function of the ribbon length L in the presence of SRI 
for several different Fermi energies E. In this case, the 
conductance decays exponentially even for single-channel 
transport. This result is similar to that previously ob- 
tained in zigzag ribbons. However, the rate of decay in 
the low-energy single-channel regime (E — 0.1 and 0.2) 
is slower than that for multi-channel transport regime 
(E > 0.4) in this case. Similar results are obtained in 
Ref but in which only short-ranged disorder at edge 
of ribbon is considered. 



B. T-matrix analysis 

The absence of localization in the single-channel region 
can be understood from the Dirac equation including the 
impurity potential term L^ mp with armchair edge bound- 
ary. To consider the amplitude of backward scattering, 
we introduce the T-matrix defined as 



T = U\ 



imp 



Us 



1 



-u h 



imp - ^ imp 

E-Ho 



(8) 



III. ELECTRONIC TRANSPORT PROPERTIES 
OF DISORDERED NANORIBBONS 



According to Ref.—, U- lmp is written as 



A. Numerical simulation 

Now we turn to the discussion of the electronic trans- 
port properties of disordered nanoribbons. We evaluate 



C'imp 



/ u A (r) u' A (r) 

u B (r) -u' B (r) 

u' A (r)* u A (r) 

V -u' B {r)* u B (r) 



(9) 
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with 



ux(r)=J29(r-Rx)ux(Rx) , (10) 
Rx 

u , x(r) = Y,9(r-Rx)c- i2K - Rx u x (Rx) , (11) 



Rx 



where ux{Rx) is the local potential due to impurities 
for X = A or B. Here g(R) with the normalization 
condition of ^2 R g{R) = 1 is the real function which 
has an appreciable amplitude in the region where \R\ 
is smaller than a few times of the lattice constant, and 
decays rapidly with increasing \R\. If only the LRI are 
present, we can approximate UA(r) = ub{t) = u(r) and 
u'a(i~) = u'g(r) = u'(r). In the case of carbon nan- 
otubes and zigzag nanoribbons, u' x(r) vanishes after the 
summation over Rx in Eq. pip since the phase factor 
e -i2K R x gtrongly oscillates in the x-direction. How- 
ever, this cancellation is not complete in an armchair 
nanoribbon because the averaging over the x-direction is 
restricted to the finite width of W . This means that we 
cannot neglect the contribution from scatterers particu- 
larly in the vicinity of the edges to u'x(r). Although 
v! x{t) becomes small after the summation, the symme- 
try of system changes for u'x(f) 7^ as we will see in the 
next subsection. 

Now we evaluate the matrix elements of £/; mp for the 
eigenstate \n, k, s) with the eigen energy of Eq. |||5J) which 
can be written as 



\n,k,s) = 



1 



\ 



-W(n,k) 



-W(n,k) 



jky 



(12) 



with the phase factor 



-W{n,k) 



ik 



y/< + fc2 



(13) 



Here it should be noted that the phase structure in 
Eq. (|12p is different between and K _ states, and 
this internal phase structures are critical for the scat- 
tering matrix elements of armchair nanoribbons as we 
discuss in the following. Using the above expression, we 
can obtain the scattering matrix element 

(n,M#imp|n / ,*V) = (ss' + e 1 ^)-^'.*'))) 

x V(n,k;n',k') , (14) 



with 



V(n,k;n',k') 



1 



, dx / dye- iik - k ' )y 

u(r) (e-^-^'^+cc.) 
(u'(r)e- i(K " +K "' )K +c.c.)] . (15) 



It should be emphasized that Eq. (| 14[) has the same form 
as that obtained for carbon nanotubes without inter- 
valley scattering (u' x (r) = 0) 18 . Interestingly, in spite 
of the fact that armchair nanoribbons inevitably suffer 
from the inter- valley scattering due to the armchair edges 
(u' x (r) 7^ 0), we can express the matrix element for the 
backward scattering as Eq. (114|) by including u' x {r) into 
V(n,k;n',k') in Eq. (| 1 5[) . This is due to the different 
phase structure between K+ and K_ in Eq. (|T2|) . 

We focus on the single-channel regime where only the 
lowest subband with n — crosses the Fermi level. From 
Eq. (Tnj, the scattering amplitude from the propagating 
state |0, k, s) to its backward state |0, —k, s) in the single- 
channel mode becomes identically zero, i.e. 



(0,-k,s\U imp \0,k,s) = 0. 



(16) 



Thus, since the lowest backward scattering matrix ele- 
ment of T-matrix vanishes, the decay of (g) in the single 
channel energy regime is extremely slow as a function of 
the ribbon length as we have seen in Fig.O However, the 
back-scattering amplitude in the second and much higher 
order does not vanish. Hence the single-channel conduc- 
tion is not exactly perfect like carbon nanotubes—, but 
nearly perfect in armchair nanoribbons. 

The present results of nearly perfect single-channel 
transport might be similar to those obtained in carbon 
nanotubes by solving the Boltzmann transport equation, 
which is valid for incoherent systems in the absence of 
inter-valley scattering—. However, our results are for the 
coherent system with inter- valley scattering by armchair 
edge and their physical mechanism is different. 



C. Symmetry consideration 

Now we give a symmetry consideration to disordered 
graphene and graphcne nanoribbons. If the inter- valley 
scattering is absent, i.e. u' x (r) = 0, the Hamiltonian 
Hq + Ui mp becomes invariant under the transformation 
of S = — i (a y ® t°J C, where C is the complex-conjugate 
operator. This operation corresponds to the special time- 
reversal operation for pseudospins within each valley, and 
supports that the system has the symplectic symmetry. 
However, in the presence of inter- valley scattering due to 
SRI, the invariance under S is broken. In this case, the 
time reversal symmetry across two valleys described by 
the operator T = (a z ® t x )C becomes relevant, which 
indicates orthogonal universality class. Thus as noted in 
Refi^S, graphene with LRI belongs to symplectic symme- 
try, but that with SRI belongs to orthogonal symmetry. 

However, in the disordered armchair ribbons, the spe- 
cial time-reversal symmetry within each valley is bro- 
ken even in the case of LRI. This is because u' x (r) 7^ 
as we have seen in the previous subsection. Thus, ir- 
respective of the range of impurities, the armchair rib- 
bons are classified into orthogonal universality class. Ac- 
tually, the application of magnetic field shows rather 
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strong magnetic field dependence on the inverse local- 
ization length in the regime of the weak magnetic field 
for both LRI and SRI cases (Fig. [3]). This is consistent 
with the behavior of orthogonal universality class. Here, 
the inverse localization length is evaluated by identify- 
ing exp(ln g) = exp(— L/£). Since the disordered zigzag 
nanoribbons are classified into unitary class for LRI but 
orthogonal class for SRIi2, it should be noted that the 
universality crossover in nanographene system can occur 
not only due to the range of impurities but also due to 
the edge boundary conditions. 

IV. SUMMARY 

In this work, we have numerically investigated the elec- 
tronic transport in disordered armchair nanoribbons in 
the presence of short- and long-ranged impurities. In 
spite of the lack of well-separated two valley structures, 
the single-channel transport subjected to long-ranged im- 
purities shows nearly perfect transmission, where the 



backward scattering matrix elements in the lowest or- 
der vanish as a manifestation of internal phase struc- 
tures of the wavefunction. These results are contrast with 
the mechanism of perfectly conducting channel in disor- 
dered zigzag nanoribbons and metallic nanotubes where 
the well separation between two non-equivalent Dirac 
points is essential to suppress the inter-valley scattering. 
Within the Born approximation, this cancellation is sat- 
isfied even for SRI, which can be clearly seen in Fig. [2Kb) . 
The dependency of conductance on the Fermi energy con- 
firms our calculation about the difference between single- 
and multi-channel transport properties. The symmetry 
consideration classifies the armchair nanoribbons into or- 
thogonal class. 
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